Dynamics of nonlinear resonances in Hamiltonian systems 
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It is well known that the dynamics of a Hamiltonian system depends crucially on whether or not 
it possesses nonlinear resonances. In the generic case, the set of nonlinear resonances consists of 
independent clusters of resonantly interacting modes, described by a few low-dimensional dynamical 
systems. We show that 1) most frequently met clusters are described by integrable dynamical sys- 
tems, and 2) construction of clusters can be used as the base for the Clipping method, substantially 
more effective for these systems than the Galerkin method. The results can be used directly for 
system with cubic Hamiltonian. 

PACS numbers: 47.10.Df, 47.10.Fg, 02.70.Dh 



1. Introduction. A notion of resonance runs through 
all our life. Without resonance we wouldn't have radio, 
television, music, etc. The general properties of linear 
resonances are quite well-known; their nonlinear counter- 
part is substantially less studied though interest in un- 
derstanding nonlinear resonances is enormous. Famous 
experiments of Tesla show how disastrous resonances can 
be: he studied experimentally vibrations of an iron col- 
umn which ran downward into the foundation of the 
building, and caused sort of a small earthquake in Man- 
hattan, with smashed windows and swayed buildings 
Another example is Tacoma Narrows Bridge which tore 
itself apart and collapsed (in 1940) under a wind of only 
42 mph, though designed for winds of 120 mph. 

Nonlinear resonances are ubiquitous in physics. Eu- 
ler equations, regarded with various boundary conditions 
and specific values of some parameters, describe an enor- 
mous number of nonlinear dispersive wave systems (cap- 
illary waves, surface water waves, atmospheric planetary 
waves, drift waves in plasma, etc.) all possessing non- 
linear resonances 0- Nonlinear resonances appear in a 
great amount of typical mechanical systems such as an 
infinite straight bar, a circular ring, and a flat plate 
The so-called "nonlinear resonance jump" , important for 
the analysis of a turbine governor positioning system of 
hydroelectric power plants, can cause severe damage to 
the mechanical, hydraulic and electrical systems It 
was recently established that nonlinear resonance is the 
dominant mechanism behind outer ionization and energy 
absorption in near infrared laser-driven rare-gas or metal 
clusters The characteristic resonant frequencies ob- 
served in accretion disks allow astronomers to determine 
whether the object is a black hole, a neutron star, or a 
quark star Q. Thermally induced variations of the he- 
lium dielectric permittivity in superconductors are due to 
microwave nonlinear resonances [7| . Temporal processing 
in the central auditory nervous system analyzes sounds 
using networks of nonlinear neural resonators 0]. The 
non-linear resonant response of biological tissue to the 
action of an electromagnetic field is used to investigate 
cases of suspected disease or cancer 0. 



The very special role of resonant solutions of nonlin- 
ear ordinary differential equations (ODEs) has been first 
investigated by Poincare at the end of the 19th century. 
Poincare proved that if a nonlinear ODE has no reso- 
nance solutions, then it can be linearized by an invertible 
change of variables (for details see [l3| and refs. therein). 
This simplifies both analytical and numerical investiga- 
tions of the original nonlinear equation, allows for the 
introduction of corresponding normal forms of ODEs, 
etc. In the middle of the 20th century, Poincare's ap- 
proach has been generalized to the case of nonlinear par- 
tial differential equations (PDEs) yielding what is nowa- 
days known as KAM theory ([lH-[lE] and others). This 
theory allows us to transform a nonlinear dispersive PDE 
into a Hamiltonian equation of motion in Fourier space 



(1) 



where Ok is the amplitude of the Fourier mode corre- 
sponding to the wavevector k and the Hamiltonian Ti. 
is represented as an expansion in powers Ti.j which are 
proportional to the product of j amplitudes a^. In this 
Letter we are going to consider expansions of Hamilto- 
nians up to third order in wave amplitude, i.e. a cubic 
Hamiltonian of the form 



n. = 



^230^iO2a3<523 + complex conj. 



ki,k2,k3 

where for brevity we introduced the notation aj = 
and = (5(ki — k2 — ka) is the Kronecker symbol. If 
7^3 7^ 0, three- wave process is dominant and the main 
contribution to the nonlinear evolution comes from the 
waves satisfying the following resonance conditions: 



I w(ki) +uj{k2) -ujiks) 

I ki + k2 - k3 = 0, 



where tt'(k) is a dispersion relation for the linear wave 
frequency and > is called resonance width. 
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If > 0, the equation of motion ([T|) turns into 

ki,k2 

If = 0, the equation of motion ([T|) turns into 

«-Sk = Eki,k2 (V'i^BiB2'5S'2'5(t^k - c^i - UJ2) 

+2V^^*BiB^6l^S{LJi - cok - LJ2)). (3) 

The co-existence of these two substantially different 
types of wave interactions, described by Eqs.([2]) and 

has been observed in numerical simulations 
and proven analytically in the frame of the kinematic 
two-layer model of laminated turbulence [isj . Dynamics 
of the layer ^ is described by wave kinetic equations 
and is well studied [l^. Dynamics of the layer ^ 
is practically not studied though a lot of preliminary 
results is already known. Namely, the layer ([3]) is 
described by a few independent wave clusters forme d by 
the waves which are in exact nonlinear resonance [l9|. 
The corresponding solutions of ^ can be computed by 
the specially developed q-class method, presented in [20| 
and implemented in [2lj . The general form of dynamical 
systems describing resonant clusters can also be found 
algorithmically [24I, as well as coefficients of dynamical 
systems [2^. Moreover, as it was demonstrated in [2^ 
(numerically) and in [25| (analytically), these clusters 
"survive" for small enough but non-zero fl. The main 
goal of this Letter is to study the dynamics of the most 
frequently met resonant clusters. 

2. Clusters. In this Letter we present some analyti- 
cal and numerical results for the three most commonly 
met dynamical systems corresponding to non-isomorphic 
clusters of nonlinear resonances - a triad, a kite, and a 
butterfly consisting of 3, 4 and 5 complex variables cor- 
respondingly. 

The dynamical system for a triad has the form 

Bi = Z B2B^, B2 = Z B^B^, B^ = — ZBiB2- (4) 

A kite consists of two triads a and b, with wave ampli- 
tudes Bja, Bjb, j = Ii2j3, connected via two common 
modes. Analogously to |26| . one can point out 4 types 
of kites according to the properties of connecting modes. 
For our considerations, this is not important: the general 
method to study integrability of kites will be the same. 
For the concreteness of presentation, in this Letter a kite 
with Bia = -B16 and ,620 = ^26 has been chosen: 

Bla = B2a{ZaB^a + Z^Bsb) , 
< B2a = Bl^{ZaB3a + Z^B^b) , (5) 
^B^a = —ZaBiaB2a , B^b = —ZhBiaB2a ■ 

A butterfly consists of two triads a and b, with wave 
amplitudes Bja, Bjt, j = 1,2,3, connected via one 



common mode. As it was shown in [26[, there exist 3 
different types of butterflies, according to the choice of 
the connecting mode. Let us take, for instance, Bia = 
_Bif,(= Bi). The corresponding dynamical system is then 
as follows: 

Bl = ZaB^aBsa + ZbB^fjB^h , 

B2a = ZaBlB^a , B2b = Z^BlB^b , (6) 

Bza = —ZaBiB2a , B^h = —ZhBiB2b ■ 

3. Integrability of resonance clusters. From here 
on, general notations and terminology will follow Olver's 
book [23] . We use hereafter Einstein convention on re- 
peated indices and f,i = df /dx' . Consider a general iV- 
dimensional system of autonomous evolution equations 
of the form: 

^{t) = ^\x^{t)), ^ = l,...,N. (7) 

Any scalar function f{x^,t) that satisfies 
j-^ {f{x\t),t)) = f / + AV,, = is called a con- 
servation law in [27t . It is easy to see that this definition 
gives us two types of conservation laws. The first type 
is the standard notion used in classical physics: the 
conservation law is of the form f{x^), i.e. it does not 
depend explicitly on time. The second type looks more 
like a mathematical trick: it is of the form / (a;*, t), where 
the time dependence is explicit. In this Letter we will be 
interested in both types of conservation law. We claim 
that they are both physically important and to show 
that we present an illustrative example from Mechanics. 
Let us regard the damped harmonic oscillator. The 
equations of motion in non-dimensional form can be 
written as: 

q=p^ p = -q-ap, (8) 

where a > is the damping coefficient. If this coeffi- 
cient is equal to zero, a = 0, then the total energy of the 
system 

Eiq,p)^ 1/2 (p^ + q^) 

is a conservation law of the first type. If a > 0, then the 
system does not conserve the energy anymore but one can 
still define a conserved quantity which is a generalization 
of E for the case a > 0. This new quantity 

F[q,p,t) = exp(Q!t) (p^ + + apq) 

is a conservation law of the second type. This means that 
for an arbitrary solution q{t),p{t) of the system dH) we 
have 

[exp(ai) {p{tf + q{tf + a Pit) g(t))] = 
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These two types of conservation law are very different 
but complementary. While the first type, E{q,p) in the 
case a = 0, defines where the motion takes place, the 
second type, F(q,p,t) in the case a > 0, defines how the 
motion takes place. In other words, the first type defines 
orbits of the dynamical system ^ and the second type 
defines its motion within the orbit. 

To keep in mind the difference between these two types 
of conservation laws, we will call the first type just a 
conservation law (CL), and we will call the second type 
a dynamical invariant. 

We say that Sys.([2l) is integrable if there are N 
functionally independent dynamical invariants. Ob- 
viously, if Sys.Q possesses {N — 1) functionally 
independent CLs, then it is constrained to move along 
a 1-dimensional manifold, and the way it moves is 
dictated by 1 dynamical invariant. This dynamical 
invariant can be obtained from the knowledge of the 
(TV — 1) CLs and the explicit form of the Sys.l]?]), i.e. 
Sys.(l7|) is integrable then. It follows from the Theo- 
rem below that in many cases the knowledge of only 
{N ~2) CLs is enough for the integrability of the Sys.([7]). 

Theorem. Let us assume that the system ^ pos- 
sesses a standard Liouville volume density 

and [N — 2) functionally independent CLs, 
H^,...,H^^^. Then a new CL can be constructed, 
which is functionally independent of the original ones, 
and therefore the system is integrable. 

The (lengthy) proof follows from the existence of a 
Poisson bracket for the original Sys.Q and is an 
extension of the general approach used in (2^ for three 
dimensional first order autonomous equations. The 
proof is constructive and allows us to find the explicit 
form of a new CL for dynamical systems of the form 
dH),®,®, etc. 

In the examples below, we always need to eliminate so- 
called slave phases, which corresponds to the well-known 
order reduction in Hamiltonian systems [2^. The num- 
ber N used below corresponds to the effective number 
of degrees of freedom after this reduction has been per- 
formed. 

Integrability of a triad, dynamical system (|4|) , is a well- 
known fact (e.g. [sot ) ^^'^ conservation laws are 

l23^W + \B3\^ /i3 = |BiP + |i33p. 

Sys.® has been used for a preliminary check of our 
method; in this case N = A. The method can thus be 
applied and we obtain the following CL: 

Lt - Im(BiB2S3*) , 



together with the time-dependent dynamical invariant of 
the form: 




Here F is the elliptic integral of the first kind, i?i < 
i?2 < ^3 are the three real roots of the polynomial 

= 2/27- (27/|- 

(/l3 + /23 ) (/l3 - 2/23 ) (/23 - 2/i3 ) ) /27(/23 - 1,^1^, + ll:,f'^ 

and 

V = - (2/i3 - J23 + {II^ - /13/23 + I^Y'^)/^ 

is always within the interval [i?2, ^3] 9 0. 

A kite, dynamical system (O, is also an integrable sys- 
tem. Indeed, after reduction of slave variables the sys- 
tem corresponds to iV = 6 and has 5 CLs (2 linear, 2 
quadratic, 1 cubic): 

{Lr = Rc(Zbi?3a — ZaB^b), Lj = lm{ZhB3a — ZaB^b) , 
hab^\Bia\' + \B3a\^ + \B3h\\ 
l2ab^\B2a\' + \B3a\^ + \B3h\\ 
Ik = lm{BiaB2a{ZaBl^ + ZbBl^)) , 

with a dynamical invariant that is essentially the same 
as for a triad, Sq, after replacing Z = Za + Z),, It = 
Ik{ZI + ZI)/Z\ /i3 = hab{Zl + Zl)/Z-'-[Ll + L])/Z\ 

l23 - l2at{Z^a+ Z^) / Z' - {I^l +L])/Z\ 

The dynamics of a butterfly is governed by Eqs.® and 
its 4 CLs (3 quadratic and 1 cubic) can easily be obtained: 

hza ~ \B2a\'^ + |-B3aP , hzb = |-B2hP -|- |i?36p j 

Iab^\Bi\' + \B:ia? + \B^,\\ (9) 

/q = ln\{ZaBiB2aBla + ZhBiB2bB^^) , 

while a Liouville volume density is p = 1. Notice that all 
cubic CLs are canonical Hamiltonians for the respective 
triad, kite and butterfly systems. From now on we con- 
sider the butterfly case when no amplitude is identically 
zero; otherwise the system would become integrable. 

The use of standard amplitude-phase representation 
Bj — exp(i0j) of the complex amplitudes Bj in terms 
of real amplitudes Cj and phases 9j shows immediately 
that only two phase combinations are important: 

fa = Sia + 02a — ^3a, fb = dlb + (^2b — ^36, 

a- and 6-triad phases (with the requirement 9ia = On, 
which correspond to the chosen resonance condition). 
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FIG. 1: Color online. To facilitate view, color hue of the plot is a linear function of time t, varying from to 1 as t runs through 
one period. Left panel: Za ^ Zt = 10/W0{integrablecase), lab ~ 2.1608. Middle panel; Za = 8/100,^6 = 12/100, /at ~ 
2.2088. Right panel: Z^ = 9/100, Zt = 11/100, hb ~ 2.1846. 



This reduces five complex equations ^ to only four real 
ones: 



dt 

dCsb 
dt 

difia 

dt 

dipt 



a COS ; 

-2'b(7iC2hCos ipb , 



(10) 
(11) 



sm ipi, 



(13) 



dt \ Czb C2b ) 

The cubic CL reads 

/o = Ci [ZaC2aCza Sm Lpa + ZbC2bCzh Sm iph) (14) 

in terms of the amplitudes and phases. This means that 
the dynamics of a butterfly cluster is, in the generic case, 
confined to a 3-dimensional manifold. Below we regard 
a few particular cases in which Sys.(I5]) is integrable. 
Example 1: Real amplitudes, = (/St = 0. In this 
case the Hamiltonian /q becomes identically zero while 
Liouville density in coordinates Csa, Csf, is p(C3q, Csh) = 
1/(71(720^26. So in this case the equations for the un- 
known CL H{C3a,C3b) are: 



pA 



C3a 



Za d 



-H, pA^ 



Zb 



d 



C2b dCzb C2a dC^a 

and from Eqs.Q we readily obtain 

H{C3a,C3b) = Zh arctan {C3a/C2a)~Za arctan {C3b/C2b) , 

i.e. Sys.® is integrable in this case. Of course, this case 
is degenerate for Iq = yields no constraint on the re- 
maining independent variables C^a , C^b satisfying equa- 
tions TO . (fTT|l. 

General change of coordinates. Going back to 
Eas. (fTUl) - (|14|) . we can jump from this degenerate case to 
a more generic case by defining new coordinates which 



are suggested by H{C3a,C3b)- The new coordinates are 
to replace the amplitudes C^a, C^b- 

aa = arctan ((73a/(72a) , ab = arctan ((73fc/(72h) . 

We choose the inverse transformation to be 

J C2a = Vhsa C0s(Q;a) , = \J hza sin(aa) 



(726 = \//23bCOs(ab) , (736 = \/ hzb sin(Q;6) , 



(15) 



so that the domain for the new variables is < eta < 
7r/2, < q;6 < 7r/2. For these new coordinates, the 
evolution equations simplify enormously: 



= -ZaCiC0S</P„ ^ = -Z6CiCOS¥P6, 

^ = ZaC\ (cotaa - tan Q!a) sincpa 



(16) 



^ = Z6(7i (cota6 - tan a6) sin .^6 - (^fp-, 
where the amplitude (7i > is obtained using eqs. 



Ci 



I: 



23a sin aa - hab sin ab 



(17) 



and the cubic CL is now 
(7i 

/o = — {Zal23a sin(2aa) sin(i^a) + ^6-^236 sin(2a6) sin(v36)) ■ 

(18) 

Equations (fTCl) - P^ represent the final form of our 3- 
dimcnsional general system. 

Example 2: Complex amplitudes, Iq — 0. Here, we 
just impose the condition /q = but the phases are 
otherwise arbitrary: this case is therefore not degener- 
ate anymore and we have a 3-dimensional system which 
requires the existence of only 1 CL in order to be in- 
tegrable: a CL is Aa = sin(2ao) sin((/3a), which can be 
deduced from Eqs. (fT6|) . Making use of the Theorem, one 
can find another CL for this case: Hneu)((73a, (73b) = (1 + 



Zb/Za) arccos 



— (1 + Za/Zb) arccos 
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FIG. 2: Color online. Chain of 4 triads, real amplitudes for 
all 9 modes. Time and amplitudes in non-dimensioned units. 



Obviously Aa and Hnew are functionally independent, i.e. 
the case /q = is integrable. 

Example 3: Complex amplitudes, Za = Zb- In this 
case a new CL has the form 



Za 



-E 



'^la'^la + ClbClb + 2(720(^306*266^36 C0s((^a - (fb) 



^2 



(^ia ^2b + (^3b) J 



(19) 



which is functionally independent of the other known 
constants of motion. Therefore, according to the The- 
orem the case Zb = Za is integrable. 

The numerical scheme is programmed in Mathematica 
with stiffness- switching method in single precision. 
For arbitrary Za, Zb the scheme has been checked by 
computing /q from ea. ([T5)) at all consequent time steps; 
there is no noticeable change of Iq up to machine 
precision. 

4. Numerical simulations. To investigate the 
general behavior of a butterfly cluster with Za ^ Zb, 
we integrated directly Eqs. p^ . fTT)) with Iq computed 
from Eq. p^ evaluated at t = and used to check 
numerical scheme afterwards. Some results of the 
simulations with Svs. ([T6|) are presented in Fig[T] Initial 
conditions aa(0) = 78/100, a6(0) = 60/100, V3a(0) = 
147/100, (pb{0) = 127/100 and values of the constants 
of motion lo = 11/2000, /23a = 4/100, /236 = 4/100 are 
the same for all three parts of Fig[T] 3D parametric 
plots are shown in space {aa,^a,'Pb) with color hue 
depending on time, so the plots are effectively 4D. The 
main goal of this series of numerical simulations was 
to study changes in the dynamics of a butterfly cluster 
according to the magnitude of the ratio C = Za/Zb- In 
FiglH left panel Q = 1, middle panel = 2/3 and 
right panel Cr = 9/11. As it was shown above, the case 
= 1 is integrable, and one can see the closed trajectory 
with period T; sa 21.7. Quite unexpectedly, rational 
Cm,Cr produce what appear to be periodic motions and 
closed trajectories with periods Tm ~ 53 and « 215 
correspondingly. A few dozen of simulations made with 
different rational ratios C = Za/Zb show that periodicity 
depends - or is even defined by - the commensurability 
of the coefficients Za and Zb- FiglT] shows that = 2/3 



gives 2 spikes in one direction and 3 spikes in the 
perpendicular direction, while (^r = 9/11 gives 9 and 11 
spikes correspondingly; and so on. This does not mean, 
of course, that in the general case Za ^ Zb solutions 
of Sys.® are periodic. But the opposite can only be 
proven analytically for in numerical simulations any 
choice of the ratio ^ necessarily turns into a rational 
number. Some preliminary series of simulations have 
been performed with the chains of triads with connection 
types as in ([6]), with maximum number of triads in 
a chain being 8 (corresponds to 17 modes). For our 
numerical simulations, resonant clusters of spherical 
Rossby waves were taken, with initial (non-dimensioned) 
energies of the order of measured atmospheric data, as 



3l| . The dynamical system for computations was 



taken in the original variables Bj. The case of real 
amplitudes is shown in Fig[2] all resonant modes behave 
(almost) periodically. Non-dimensioned units for time 
and amplitudes were chosen to illustrate clearly the 
characteristic behavior of the amplitudes. 

5. Galerkin method versus Clipping method. 

Speaking very generally, Galerkin method (GM) allows 
one to reduce infinite dimensional systems, described 
by PDEs, to low dimensional systems of ODEs. Three 
steps have to be performed: a) choice of ansatz functions 
(modes); b) choice of number of modes A^; c) construc- 
tion of the resulting reduced system. While in the case of 
a linear PDE, Fourier modes is the natural choice of the 
ansatz functions, the number of modes is mostly defined 
by the available computer facilities and the construction 
of the corresponding ODE is obvious, the case of a non- 
linear PDE is much more involved |32l | . 

The Clipping method (CM) has been introduced in 
[33! in order to deal with evolutionary dispersive nonlin- 
ear PDEs. The idea of the CM is very simple: under 
some physically relevant conditions, it is enough to re- 
gard the dynamics of the resonantly interacting modes. 
These resonant modes can be found systematically us- 
ing the q-class method [2QI|. All other modes behave as 
linear and can be clipped out. Numerical evidence of 
the effectiveness of this approach was presented in [23 |. 
using a pseudo-spectral model of the barotropic vortic- 
ity equation. There, it was observed that the most en- 
ergetically active modes were the resonant modes (221 
modes were excited, among them 39 resonant). These 
resonant modes appear in clusters, i.e. a small number 
of modes, and the energy of each cluster is conserved (see 
24|, Fig.l). As for each of the non- resonant modes, their 
energies are approximately constant during many periods 
of energy exchange of the resonant modes ([13], Fig. 3). 

Implementation of the q-class method for various 
dispersion functions |2l| . together with explicit con- 
struction of the dynamical systems corresponding to 
the exact resonances [2^, allows one to eliminate the 
arbitrariness of the choices at all three steps a), b) and 
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c) of the GM. Instead of one big system of ODEs of order 
N, we have a few independent systems of ODEs, of order 
N <^ N, and these systems are often integrable. Another 
example can be found in 22]: ocean planetary motions, 
128 resonant modes among 2500 Fourier harmonics 
in the chosen spectral domain, all clusters and their 
dynamical systems are written out explicitly. In dozens 
of studied 3- wave systems, the amount of resonant modes 
does not exceed 20% of all modes in the spectral domain. 

6. Conclusions. Our analysis and general mathemati- 
cal results [l^ on resonant clusters are valid for arbitrary 
Hamiltonian Tij^j > 3, though computation of clusters 
in the case j > 3 is more involved 21 1. 

Clipping method has at least three advantages com- 
pared to Galerkin truncation: i) Numerical schemes in 
CM can be truncated at a substantially higher wavenum- 
ber than in GM, depending not on the computer facilities 
but on some physically relevant parameters (say, dissipa- 
tion range of wavenumbers). ii) Most of the resulting 
dynamical systems corresponding to each resonant clus- 
ter is analytically integrable. iii) The solutions obtained 
from the integrable cases could be used to parameterize 
the numerical solutions of non-integrable systems found 
for bigger clusters. This work is in progress. 

Last not least. Even for one specific PDE, it is a 
highly non-trivial task to prove that Galerkin truncation 
is a Hamiltonian system and to construct additional 
conserved quantity 3J| . Clipping method combined with 
the constructive procedure based on the Theorem above 
(Section 3), allows us to produce physically relevant 
dynamical systems and to find additional conservation 
laws systematically, for a wide class of evolutionary 
dispersive nonlinear PDEs. 
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